Many-body fits of phase-equivalent effective interactions 
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In many-body theory it is often useful to renormalize short-distance, high-momentum components 
of an interaction via unitary transformations. Such transformations preserve the on-shell physical 
observables of the two-body system (mostly phase-shifts, hence unitarily-connected effective inter- 
actions are often called phase-equivalent), while modifying off-shell T-matrix elements influential in 
many-body systems. In this paper I lay out a general and systematic approach for controlling the 
off-shell behavior of an effective interaction, which can be adjusted to many-body properties, and 
present an application to trapped fermions at the unitary limit. 
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PACS numbers: 21.60.Cs,21.30.Fc,67.85.Bc 



The force between two particles has to be determined 
empirically, but some matrix elements are easier to de- 
termine than others. For an isolated two-body system 
one can measure elastic phase shifts and bound state 
eigenvalues, but not off-shell (inelastic) matrix elements 
(specifically, the elements of the T-matrix). Such ma- 
trix elements show up when the two-body system is not 
isolated, i.e., embedded in a many-body system, but dis- 
entangling the off-shell two-body matrix elements from 
many-body data is not possible. 

Following modern trends in effective interaction the- 
ory [IHE], I will argue that one can flip the ambiguity of 
off-shell matrix elements into an advantage. This paper 
illustrates a general method to find the unitary transfor- 
mation that preserves two-body observables while simul- 
taneously providing the best fit to many-body data. 

For example, despite much effort we do not have 
a unique prescription for the force between nucleons. 
Nonetheless high-precision data on two-nucleon systems 
Q strongly constrain any description, leading to a vari- 
ety of competing interactions [7HlO|| which all, by con- 
struction, have indistinguishable on-shell T-matrix ele- 
ments; these are called phase-shift equivalent potentials. 
The off-shell T-matrix elements, which do differ, can only 
be compared through many-body observables, including 
binding energies and excitation spectra. 

In a parallel vein, consider configuration-interaction 
(CI) many-body calculations, which have several advan- 
tages but are vulnerable to strong short-range, high- 
momentum components of the interaction, because prac- 
tical considerations require truncation, indicated by a 
projection operator P. The problem is the eigenvalues 
of PHP are not the same as the eigenvalues of H in an 
infinite or even very large space, and converge slowly with 
increasing dimension of P, a particularly severe problem 
for CI calculations. Therefore one turns to effective in- 
teractions, generated using a unitary transformation 

H cS = U^HU. (1) 

The hope is that the eigenvalues of PH c gP in the trun- 
cated space converge to those of H in an infinite space. 
One common strategy is the so-called 'cluster approxi- 
mation': forcing the eigenvalues in the truncated space 



to agree with the 'exact' values of the original space, but 
only for two- (or on occasion three-) particle systems, and 
then apply to many-body systems. 

Modern, rigorous effective interaction methods explic- 
itly or implictly apply a unitary transform that leave 
the on-shell behavior unchanged but which dial away 
the troublesome high-momentum part of the interaction. 
These methods, which include but are not limited to 
Okubo-Lee-Suzuki (OLS) ll|, the unitary correlation op- 
erator method (UCOM) 2], the similarity renormaliza- 
tion group (SRG) Q, and Alhassid-Bertsch-Fang (ABF) 
IH, generate a portfolio of new 'phase-equivalent' inter- 
actions with the same on-shell T-matrices but different 
off-shell matrix elements. 

The two-body cluster approximation is complicated by 
the fact that three- body forces naturally arise out of effec- 
tive field theory [Hf- In practice, however, three-body 
forces make for computational difficulties. Thus several 
authors have sought to minimize three-body forces by 
exploiting the interplay between three-body forces and 
off-shell matrix elements 0, [H, EII • 

This paper combines several of these ideas. Borrowing 
from phenomenological fitting of interactions [3, HH , I 
show how one can choose a general unitary transforma- 
tion to generate a phase-equivalent potential which is a 
"best fit" to many-body data. 

To illustrate, I take on a specific challenge: spin-1/2 
fermions at the unitary limit, the so-called "Bertsch prob- 
lem," in the context of an external harmonic trap [l6l[l7| . 
The interaction is zero range and has infinite scattering 
length, similar to the short-range, large scattering length 
nuclear force. The Hamiltonian is 
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The ground state energy for the 3-body case is known 
analytically [HI, and for N > 4 I take as 'exact' the re- 
sults from correlated Gaussian and fixed-node diffusion 
Monte Carlo calculations [TtJ (with quoted statistical and 
systematic errors of only a few percent, which I leave 
out). Table I lists the energies adopted in this study. 
For purposes of comparison I group the energies into two 
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sets: Set I includes the ground states and selected ex- 
cited states for N = 3,4, while Set II is comprised of the 
ground state energies for N = 3-10. 



TABLE I: Adopted energy levels [13, [13] in units of HQ of the 
trapping potential and assignment to comparison sets. 



N 


V ; S energy 


Set 


3 


\-;\ 4.27 
0+; | 4.66 


1,11 
I 


4 


+ ;0 5.05 
2+;0 5.91 
1 + ;1 6.58 


1,11 
I 
I 


5 


l";i 7.53 


II 


6 


+ ;0 8.48 


II 


7 


l + ;± 11.36 


II 


8 


0+;0 12.58 


II 


9 


+ ;± 15.69 


II 


10 


+ ;0 16.80 


II 



The Hamiltonian © can be separated into the center 
of mass plus the relative Hamiltonian: 

ffrel = ~ + \^ 2 r 2 - V S (r) (3) 

where f = r\ — ri is the relative coordinate and fx = m/2 
the reduced mass. For application to configuration in- 
teraction (CI) calculations, one computes the matrix ele- 
ments of H Te \ in a harmonic oscillator basis, (n' l\H TC i\nl) 
(which are only nontrivial for the relative s-wave or 
1 = 0), and then transform to the lab frame two-body 
matrix elements via Brody-Moshinsky brackets [l9| . 

I truncated the relative space to n cu toff = 5, which 
means for (n'l\H ie i\nl) I use n,n' — 1, . . . , 5; this corre- 
sponds to including up 8h£l in excitation energy in the 
relative space. (Calculations with different n cu toff had 
similar results.) For the lab frame, also using a harmonic 
oscillator basis, I truncated the single-particle space to 
four major shells, that is, 0s, Op, IsOd, and lpO/. 

The ^-interaction must be regularized to fix the scat- 
tering length; for CI calculations in an oscillator basis, 
the interaction strength Vq depends on n cu t ff, the num- 
ber of s-wave basis states used . For the 'bare' interac- 
tion, fixing the scattering length is equivalent to fixing 
the grou nd state energy of the relative two-body state 
[H SH3| ; due to truncation, the excitation energies differ 
from the infinite space values. Alternately, I also used, 
in the n cuto ff = 5 space, the ABF interaction, which fixes 
all five eigenenergies of H re \ in the truncated space to the 
correct (infinite space) values. 

Going to three or more particles immediately illus- 
trates some of the headaches of effective interactions, as 
can be seen in Table II. Consider the rms error between 
calculated and target (experimental) energies, 
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For 'bare' interactions, the rms error on Set II (ground 
state energies, in units of Ml of the trapping potential, for 
A=Z-10) is 1.16, while for ABF it is 2.32. That is, forcing 
the effective interaction to have the correct eigenvalues 
for the two-body system can lead to larger errors in the 
many-body system. (This comparison, oversimplifies the 
story; please read the original [J].) 



TABLE II: Root-mean-square error between adopted exact 
energies (Table I) and calculated CI energies. 'Starting H le \ 
refers to using either the 'bare' or ABF regularized interaction 
(see text). 'Generators' refer to the set of operators used in 
the unitary transformation to minimize the rms error: 'none' 
means no transformation was performed, 'all' means all ten 
generators were used, and d/dr denotes using only the single 
generator d/dr in the relative space. Units are HQ of trapping 
potential. The unitary transformations were fit to either Set 
I or II; '/ — > //' refers to fitting to Set I, but computing the 
rms error on Set II. 
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The apparent paradox of a bare interaction doing bet- 
ter than a renormalized one can be understood by writing 
U = exp(A), where A is an anti-Hermitian two-body op- 
erator. Then U induces many-body terms in H c g cutoff 
by the two-body cluster approximation. 

Previous authors made specific implicit choices in the 
form of their off-shell behavior to minimize the effect of 
three-body interactions. The strength of the correlation 
operator in UCOM [2j, which shifts nucleons away from 
each other, is fit to minimize the error in the A = 3 and 
4 ground state energies. The inside-nonlocal, outside- 
Yukawa (INOY) potential [l2j is adjusted to both two- 
body data plus the triton binding energy, while the J- 
matrix inverse scattering potential (JISP16) is fit to two- 
body data plus binding energies up to 16 O (T3 |. 

I propose a general method to choose the 'best' unitary 
transformation. Expand A in a set of generators 
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The dependence of the calculated many-body eigenener- 
gies {E a (c)} on the Ci can be expanded to first order 



E a (c]^E a {0) + J2^Ci 



(6) 



The required derivatives can be found via the Hellman- 
Feynman theorem pll ]: 
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da 
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FIG. 1: (Color online) Comparison of ground state energy 
for N particles in a harmonic trap, interacting via a contact 
interaction: (black) circles are "bare" results (fixed in rela- 
tive frame at n cuto g = 5), (red) squares are the best fit of all 
generators of a unitary transformation, and (blue) triangles 
are the exact results. All energies in units of ftSl of the trap- 
ping potential. Slight horizontal offshifts were introduce to 
aid visibility. 



Now consider a target set of (experimental) many-body 
energies {E^}- By minimizing the root- mean-square dif- 
ference (0]) between calculated and target (experimental) 
energies, one obtains 
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dE a dE Q 
dci dcj 



dE a 
da 



(8) 



which can be solved for the Cj . Because Eq. is only a 
linear approximation, one may need to iterate to converge 
on a best solution. This methodology is very similar to 
that used to fit two-body matrix elements to low energy 
nuclear spectra [H], [H| , except that in those cases one di- 
rectly fits Hamiltonian matrix elements in the lab frame, 
and here I fit the coefficients of generators of a unitary 
transformation (and in my application below I work in 
the relative frame). 

(As an important technical point, the matrix My = 
J2 a i^ri^r 1S 111 g enera l nearly singular, so one uses 
singular value decomposition [22l [23| to find the domi- 
nant modes. Here singular value decomposition is noth- 
ing more than using a spectral (eigenvalue) representa- 
tion of My and keeping only nonsingular terms. ) 

This methodology is flexible and general. One can use 
as many or as few generators as desired or thought phys- 
ically relevant, and can constrain to a best fit using an 
arbitrary choice of many-body energies. 

I applied this prescription to my example system. 
As described above, in the relative space (n',l — 
0\H IC \\n, 1 = 0) with n cuto ff = 5 is a 5 x 5 real symmetric 
matrix (using either the 'bare' or ABF matrix elements). 
I then introduced a general real orthogonal transforma- 
tion as in Eq. (JTJ, where U also is a 5 x 5 matrix , with 10 
possible antisymmetric generators Aj. (I only performed 



the transformation in the relative s channel. In principle 
an effective interaction could generate nontrivial matrix 
elements in other channels. I leave this to future inves- 
tigations.) I then fit the parameters Cj in Eq. ([5]) by 
minimizing the rms error on either Set I or Set II. The 
results are shown in Table II (the energies are all in units 
of hfl of the trapping potential) and illustrated in Figure 
1. There is a dramatic reduction in the rms error. 

The best fit parameters for Sets I and II differed, of 
course. I also considered the error in extrapolation, by 
fitting to Set I but calculating the rms error on on the 
larger Set II; found in Table II in the rows marked 'I— >II,' 
the error in extrapolation is at worse fifty percent larger 
than the best fit to Set II. (The first number in these 
rows, in parentheses, is the original rms error, as there is 
no fit.) 

Interestingly, although ABF starts with a larger rms 
error, it yields fits with smaller errors. 

By insisting that the on-shell two-body matrix ele- 
ments (relative-space eigenvalues) remain invariant, and 
keeping only the s-wave relative channel, the above is the 
best fit that one can obtain. The remaining residual er- 
ror is due to induced three-body forces that cannot be 
fully replicated using a two-body force [24j . 

The SVD of My was always dominated by one or two 
eigenvalues and thus one or two generators. In addition, 
one would appreciate simple physical insights into the 
unitary transformation. I therefore tried the manifestly 
antisymmetric generator d/dr (where r here is the rela- 
tive coordinate), because of simplicity but also because 
of similarity to UCOM 0. Table II shows this single 
generator significantly reduces the rms error. 

The minimization will be model-space dependent, but 
one can see this as either a weakness or a strength. Most 
many modern interactions already have an intrinsic cutoff 
dependence. Here I have given how to best "tune" an ef- 
fective interaction to a model space. One can also choose 
key energy levels to fit to, not only ground state ener- 
gies but excited states that contain important physics 
such as spin-orbit splitting. Also, critically, it will be im- 
portant to constrain the unitary transformation to other 
observables such as the rms radius. Given that only a few 
degrees of freedom were used to minimize the energies, 
this is not impossible and is under investigation. 

To summarize, I have discussed general unitary trans- 
formations which produce effective interactions, and have 
shown how one can use many-body data to improve the 
interaction while preserving eigenvalues in the relative 
space (on-shell T- matrix elements). This methodology 
is a generalization of modern effective interaction theory 
and of previous specific attempts to reduce the need for 
three-body interactions. With a simply defined yet nu- 
merically challenging case of trapped fermions at the uni- 
tary limit, I demonstrated improvement in ground state 
energies. 

One of the claims of current methods is they provide 
reliable error estimates for predictions. While I have not 
addressed this important issue, one should see that the 
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error estimates are dependent on the choice of unitary 
transformation, an issue that has not yet been addressed 
in a deep way. At least I hope to provokes a closer inves- 
tigation of competing effective interaction methodologies 
(and their claimed error estimates): OLS, SRG, UCOM, 
and other specific choices. 

Extending this work to nuclei and to the inclusion of 



observables is under way. 
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